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The Casimir effect as a candidate of dark energy 

Jiro Matsumoto 
Department of Physics, Nagoya University, Nagoya ^64-8602, Japan 

It is known that the simply evaluated value of the zero point energy of quantum fields is 
extremely deviated from the observed value of dark energy density. In this paper, we consider 
whether the Casimir energy, which is the zero point energy brought from boundary condi- 
tions, can cause the accelerated expansion of the Universe by using proper renormalization 
method and introducing the fermions of finite temperature living in 3-|-n-|-l space-time. We 
show that the zero temperature Casimir energy and the finite temperature Casimir energy 
can explain dark energy and dark matter, respectively. 

§1. Introduction 



The discovery of the accelerating expansion of the Universe in late 1990'P has 
brought many models of dark energy in this ten and several years. Because the 
modifications of the Einstein equation are necessary to realize the accelerating ex- 
pansion when we assume that the Universe is isotropic and homogeneous. Including 
the cosmological constant A into the Einstein equation is generally regarded as the 
most possible model to explain the accelerated expansion, because of its simplicity. 
This model, however, has a so-called fine-tuning problem. There are several ways 
to illustrate this problem, the fine tuning problem means that the energy scale of 
the cosmological constant A to explain the acceleration of the Universe is too small 
compared to the Planck scale Mpi, which is the energy scale of the gravity. There 
^ , have been attempts to explain the accelerating expansion of the Universe not by a 

t I cosmological constant but by dynamical fields without such a fine tuning. And there 

' are also modified gravity models to explain such an acceleration. As a matter of 

, fact, fine tunings of the free parameters are usually included in these models. 

I On the other hand, there are attempts to explain the observed value of the cos- 

■ mological constant by the zero point energy of the quantum fields.'^' However, there 

would not have been persuasive explanations so far. In this paper, the zero point 
energy generated from boundary conditions, which is called the Casimir energy, is 
considered. To be concrete, the Casimir energy of the graviton and the fermions are 
^ \ considered when space-time has compact extra dimensions. Similar investigations 

are seen in the literature written by B. Greene and J. Levin.® The main differ- 
ence of this paper from B. Greene and J. Levin's is to introduce finite temperature 
fermions instead of fermions and scalars at zero temperature. It will be shown later 
that reasonable models of cosmology cannot be constructed only from zero temper- 
ature fermions with periodic boundary conditions and the graviton, and the finite 
temperature effect can play a role of massive scalar fields. 

In Sec. II, the finite temperature Casimir effect on E? x {S^)"^ space is calculated. 
And the cosmological influence from the Casimir effect is investigated in Sec. Ill by 
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using the result from Sec. II. In Sec. IV, the influences to the experiments when 
we assume the fermions as neutrinos are shortly discussed. Concluding remarks are 
given in Sec. V. Units of /cb = c = = 1 are used and the gravitational constant SttG 
is denoted by = Svr/Mpi^ with the Planck mass of Mpi = G^^/^ = 1.2 x lO^^GeV. 

§2. The finite temperature Casimir effect on X (S^)^ space 

In this section, we will calculate the finite temperature Casimir effect on x 
(S^)" space. We use the formulation and results in the literature written by S. Bel- 
lucci and A. A. Saharian.® The following calculations are the extensions of a part 
of the results in ilP to use them in the case of finite temperature. And formulas 
written in Integrals and Series^^ will be used in the following without provisos. 

The energy-momentum tensor of Dirac field on x space is given as follows: 



where is defined as 



^k,l„ 



\ 



i=4 



In Eq. (|2-2p . T? does not mean the summation with respect to i and ki = Inli/d 
for i > 4. The common periodic length d are used in Eq. ()2-ip and ()2-2p because 
the equivalence of each extra dimension has been assumed here. We denote the 
real degrees of freedom of particles by A^. These Eqs. ()2Tp . ()2-2p can be shown in 
the case of scalar field. So we assume Eqs. ()2-ip . ()2-2p are also held in the case of 
the graviton. The graviton and the fermions are only considered as the degrees of 
freedom of particles in the following. 

r 4+n 1 -. 

If we assume neutrinos as the fermions then N = Ny = 2L 2 J for each gener- 
ation of neutrinos and A^ = Ng = —{n + 4){n-\-l)/2 for the graviton as the degrees 
of freedom. Here, we assume that the neutrinos are Dirac type. The negativeness of 
the degrees of freedom for the graviton comes from the difference of statistics. 

Then, we find the following relations from Eqs. ()2-ip and ()2-2p : 



Po = 3|(n + l)po--^ + d^ 

where po = (0|rO|0), po = " ELi(0|7;iO)/3 and pm = - Er=4 (0|^|0)/n, because 
the signature of the metric is assumed to be , — ). Therefore, we 

do not need to calculate po and pbO; respectively. Calculations of po will be omitted 



'The Casimir effect as a candidate of dark energy' 



3 



and we only give the expression for po because the expression of po is given in 
former researched ''^''^^ and we wih conduct similar calculations to derive the finite 
temperature correction to the Casimir effect as follows, 



po = N 



4 + TI 

2 



K 4+n 
2 



4,+n 



Here, Kn{z) is the modified Bessel function of the second kind defined by 



Kn{z) 



-ztu2 



it' 



, 2n-l 

l) — dt. 



In the case of massless particles, Eq. ()2-5p is simplified to 

/3+n 



Po = N 



4+n 
2 



(2-5) 



(2-6) 



(2-7) 



By using the Matsubara formalism and the zeta function regularization, we obtain 
the following expression of the free energyP 



AtFq = ksT ^ In ('l =F e'^'^'' 



(2-8) 
(2-9) 



where J indicates the eigenvalue of the momentum and E'q®" is the Casimir energy 
of zero temperature. In the second line of Eq. ()2-9p . the above sign of =F is applied 
to bosons and the below sign of =F is applied to fermions. The expression of E'lf^^ = 



Vd^po is given in Eq. (|2-5p so that we consider the density of AtFq: 
d^k 



1 



(Fd") d''l3 J (27r)3 



(2-10) 



(i4,i5,---,/n+3)eZ" 



where y is a volume of the three dimensional space. The energy density is, however, 
singular in the expression of Eq. ()2T0p . Therefore, we define the energy density pT 
by subtracting the zero point energy without boundary conditions from Eq. (|2-10p 
as defined in the Casimir energy of zero temperature: 

d^k 



PT 



J (27r)3 



\N\ 



(27r)"+3 
d^k 



d"/3 J (27r)3 



E 

s=l 



In =F e 
~ (±1)^ 



-pyjrn^+k^" 



E 

J/4,/5,-,«n+3)GZ" 
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I' dP-k' si3y/m?+k'^+k 
J (27r)" 

^ s=l j=0 



(2-11) 



d^+^k 
(27r)3+i 



E 



(i5+j,-,'n + 3)eZ"-J- 



27r 



(2-12) 



Eq. (j2-12p can be transformed into 



A,{s,d,(3)=4s/3d^+' 

(«5+,,---^+3)eZ"-^-i u=i 



5+2 



27r 



2 



m2 + (27r/d)2 ^^+3^^. /2v/s2/32 + w2d2 



j+5 
2 



(2-13) 



The derivation of Eq. ()2-13p is written down in Appendix. Substituting Eq. ()2-13p 
into Eq. (j2-lip . we obtain 



n-l 



PT = -MN\d^-^'Y{±iyY E E 

«=i i=o ih+j,- ,«n+3)eZ"-j-i n=i 



2 



rri' 



Vs2/32 + ^,2^2 



5+2 
2 



5+3 



2tt 



(2-14) 



Eq. (|2-14p can be rewritten in more simple form as: 



PT 



2\N\m 



4+n 



K i+n 
2 



, 4+n 

(2vr) — 



■E(±ir E 

s=l (/4,-A+3)eZ" /^2s2+^2^3+n^2 



(2-15) 



The way of this transformation is identical with that in the appendix of S. Bellucci 
and A. A. Saharian. 

4) [*)] From Eqs. 1^ and I^HM . we obtain the finite temperature 



*■* To explain the way of the transformation briefly, making Z\j(s,/3,d) from Eq. (|2-15p and 
comparing it to Eq. (|2-13p gives the equivalence of Eqs. (|2T4p and (I2-15|I . 
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Casimir energy density /Ocasimir = Po + PT in the collected form: 



K 



4+n 



{i-25soT{-ir- 



(2-16) 

where the notation ( — 1)" = 1 for bosons and (—1)" = —1 for fermions is used. 
§3. Cosmic acceleration caused by the Casimir effect 

3.1. Friedmann-Lemaitre equations 

First, we assume Friedmann-Lemaitre- Robertson- Walker (FLRW) metric of the 
flat space, 

ds^ = -dt^ + a^{t)6ijdx'dx^ + b'^{t)5<y^dx^dx'^, (3-1) 

where i,j = 1, 2, 3 and 0, "v^ = 4, 5, • • • , n + 3. And ^, = 0, 1, 2, • • • , n + 3 will be 
used. Then, we can write down Friedmann-Lemaitre (FL) equations as follows: 

3^2 + -iuHHi, + Ki{n - l)Hl = K^d''{t)p, (3-2) 

-3^2 -2H - 2nHHb - ^n{n + 1)H^ - nHb = K^d^{t)p, (3-3) 

-6H^ -3H + 3(1 - n)HHb - ^n{n - 1)hI + (1 - n)Hb = K^d''{t)pb. (3-4) 

Here, d{t) = dQb{t) is the periodic length of the compact dimensions and do is the 
current value of it. The Hubble rate of the three dimensions is defined by H(t) = 
a{t)/a{t) and the expansion rate of extra dimensions is defined by Hf,{t) = h{t)/b{t). 
p, p and pb are defined as p = —Tq, p = Tl/2, and pb = /n, respectively. By 
eliminating the terms H{t) and H'^{t) from Eq. ()3-4p . the equation of motion for the 
size of extra dimensions is derived: 

On the other hand, we obtain the equation of continuity V^T^'^ = in the following 
form, 

p + m{p + p) + nHb{p + Pb) = Q. (3-6) 

The Casimir energy density, pressure and the pressure in extra dimensions of zero 
temperature are given by Eqs. ()2-4p and ()2-5p . Substituting Eq. ()2-5p into Eq. ()2-4p 
for pq yields the relation pQ = —po, correctly. However, the total pressure from the 
Casimir effect have not been obtained. It can be seen from Eqs. (|2-12p and (|2-15p 
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that the Casimir energy density in the FLRW Universe is obtained by putting d and 
(3 into d{t) and in Eq. (j2-15p . The concrete form of the will be given later. 
From the definition of the free energy, we obtain the following expressions of the 
pressures: 

f^casimir- Pca^imir 3 g^^^y ') 

-_ 6(t) f gjOcasimir , ^/Ocasimir 9(3 \ 

Pfe, casimir — Pcasimir I 00 / ' (■^'"J 



n 



9b{t) 9(3 db{t) J ' 



where pcasimir is treated as Pca.simh[bit), (3 {ait), b{t)}] instead of /Ocasimir[d(t), ■ 
These Eqs. ()3-7p and ()3-8p are consistent with equations in ()2-4p and Eq. (|3-6p . 
respectively. 

Next question is the a{t) and b{t) dependence of the temperature. Prom the 
definition of the free energy, we have the entropy: 

5cas.mir oc -^ksE^^!! oc a^6"/?^^. (3-9) 

The adiabatic expansion of the Universe implies, 

'^'S'casimir = 0- (3' 10) 

Thus, if the term proportional to (3~^~^{t) is a dominant component of pT, then we 
have, 

3 n 

(3{t) (xa^{t)b^{t). (3-11) 

3.2. Expansion history of the Universe 

It is difficult to solve Eqs. ()3-2p - ()3-4p . but we can find the qualitative behavior 
of the solution. We can realize the following scenario if the values of parameters are 
properly selected. 

As the Universe expands, the matter density decreases, and the right hand side 
of Eq. (|3-5p becomes small and it can vanish at a certain time. Then the Universe 
evolves as to make the right hand side of Eq. ()3-5p zero. To say more exactly, the 
scale factor b{t) moves to satisfy b{t) ~ and p + 2p{, — 3p ~ 0. Because there is the 
stable solution that satisfies b{t) « and p + 2pi, — 3p « 0. It is understood from the 
behavior of the perturbation 6b{t) = b{t) — isoiution- The condition 5b{t) > makes 
p + 2piy — 3p to be less than zero and yields Hi,{t) = b{t)/b{t) < 0, and 6b{t) < 
makes p + 2pf, — 3p to be more than zero and yields Hf,{t) > 0. If the conditions of 
this stable solution are satisfied, the time evolution of the typical four dimensional 
energy density of zero temperature, b~'^{t), is faster than constant and more slowly 

12 

than a~ (t), and the evolution of the typical four dimensional energy density of finite 

temperature, /3~"~ is more slowly than a 3+" (t) and faster than a~ {t). 

These width of the time dependence are caused from the difference between b{t) = 
and p+2pf, — 3p = 0. Here, (p+2pf, — 3p) casimir ~ /S"^""" is assumed and Pmatter oc 
is considered. When n — )• c«, these conditions agree with each other and the solution 
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is determined completely. It is found from the behavior of four dimensional energy 
densities of the Casimir effect that the zero temperature part of the Casimir energy 
density act as dark energy and the finite temperature part of that act as dark matter 
in this scenario. Moreover, the negativeness of the Hubble rate of extra dimensions, 
which is caused by the trapping of the stable solution, contributes to the dark energy 
density, effectively. 

We will check whether the solution satisfies such a stability when n = 1 as an 
example. The Casimir energy density is given as: 



— y 

(^,/4)7^(0,0) 



15 1 



2 



3 



+ 



m 



(3-12) 



The first line of Eq. ()3-12p are contributions from the graviton and the second line of 
Eq. ()3-13p are contributions from massive fermions, respectively. The leading terms 
of the energy density are given by, 

15/1 1 



47r2 [ yd^ md^ rn^d^ J \j3'^ m?l3^ 

(3-13) 

Eq. ()3-13p indicates that if/3<Cl/mord^l/m then pcasimir > 0, else if /? » 1/m 
and d » 1/m then /Ocasimir < 0, in the case of > 5. At the same time, the 
balancing condition, p + 2pi, — 3p, is rewritten by, 

P + ^Pb — 'iP — ^/^casimir — ^" 771 r T 1" /'matter 

act 4 op 



47r2 \d^ 4/3V 47r2 [ \ d^ d^ md"^ m'^d^ 
+ ^-'-^ (-^2+^ + 1^^ + iJ^) } + . (3-14) 

We find that if d » /3 then p + 2ph — 3p can be zero when /Ocasimir > 0, else if 
d <C /3 then p + Ipb — 3p can be Z6rO wll6n /^casimir 

< by combining Eqs. ()3-13p and 
()3-14p . Therefore, the condition, d > f3, must be satisfied to be /Ocasimir > when 
p + 2pb — 3p = 0. Then, we find that the solution is stable because the Casimir 
components of the balancing condition, p + 2pi, — 3p, is negative in the case of 6 ^> 1, 
and p + 2pi) — 3p is positive in the case of 6 <C 1. If the finite temperature effect is 
not considered, then p + 2p{, — 3p = is not realized when pcasimir is positive so that 
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dark energy cannot explained only from zero temperature fermions with periodic 
boundary conditions and the graviton. 

We obtain the evolution equation of a{t) by eliminating the terms proportional 
to Hb{t) and H^{t) from Eqs. (|3^ -(|3^ 



^ V6"F) = -^d" {p + (n - 1> - npfe} . (3-15) 



a^ft" dt^ ' 2 + n 
Combining this equation (|3-15|) and Eq. (|3-5|) gives 

1 d 



a36" dt 



{a^'^iH -Ht)}= K^d'^ip - pb). (3-16) 



On the other hand, the temperature of neutrinos are derived from the entropy con- 
servation when we assume neutrinos as the fermions which cause the Casimir effect: 

1 n 

where T^o = 2.725K, the scale factors a and b are normalized at the time of de- 
coupling, so that fldec = ^dcc- And we have only considered the leading terms of 
the entropy, therefore there appear correction terms proportional to the powers of 
mass. We can find that Tj^q of this scenario can be larger than that of ylCDM 
model, T^o^A = 1.945K. We need to be careful to the fact that the temperature given 
in Eq. (j3-17p is not the temperature of three dimensional spaces but the effective 
temperature of 3 -|- n dimensional spaces. 

As the end of this section, an example of values of the parameters, mi, n, d and (3 
is given when neutrinos are assumed as the fermions which cause the Casimir effect. 
The reason why neutrinos are considered is that the scale of masses of neutrinos 
makes the energy scale of the Casimir energy near to the observational value of dark 
energy. Here, mi is the mass of vi and the masses of 1^2 and ^3 are determined 
from Z\m| = 7.58 x 10^^ eV^ and |^m^| = 2.35 x 10~'^eV^'^ by assuming normal 
hierarchy. The appropriate values of the parameters are mi < 0.1 meV, n = 8, 
d = 66 ;um and T^q = 30K. Then dark energy density and dark matter density can 
be explained by the Casimir energy density. Here, the temperature of the graviton 
is treated as 0.6/1.945 x T^q. However, the temperature of the graviton does not 
so much affect to the energy density, because the obtained energy density from the 
Casimir effect changes within 10% for the variation of temperature, TgQ = (0 — 
0.8)/1.945 X T^o- 

In brief, the values of the parameters are determined in the following procedure. 
First, the practically free parameters are only mi and n in this model. Because d is 
determined by Eq. ()3-5p if mi, n, and (3 are fixed, and j3 is almost determined by the 
generation of the balancing point, p + 2pb — 3p ^ 0, and the condition, /Ocasimir > 0, 
on the balancing point. It is also noted that /3 is almost constant on the balancing 
point. On the other hand, mi and n are determined to fit the energy density from the 
Casimir effect on the balancing point to the values from observations. If n < 8 then 
there are no solutions to fit the cosmological observations in the range of mi G M. 
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For example, if n = 6, mi < 0.1 meV, d = 23.5 and T^o = 80K, we obtain 20 
times larger value of the energy density compared to the observed value of the critical 
density. Else if n > 8 except for n = 9 then mi, d, and /3 must to be large compared 
to those of n = 8 to fit the values from observations. 

§4. Influences to the experiments and the observations 

In this section, we consider the influences to the experiments and the observa- 
tions when we assume neutrinos cause accelerating expansion of the Universe. First, 
this model has a continuous symmetry of translation for the three dimensional space 
and a discrete symmetry of spatial translation for extra compacted spaces, because 
the coupling to the 3 + 1 dimensional matter with neutrinos breaks a continuous 
symmetry but the geometry (S*-^)" keeps a discrete symmetry. Meanwhile, the mo- 
mentum of the extra dimension is discretized by the same reason. Therefore, there 
is the momentum conservation, entirely. The momentum conservation guarantees, 




so that there is no influence to the neutrino oscillations if we consider the zero mode. 

Next, let us consider the correction for the decay width of the Z boson. The 
effect of extra dimensions are expressed as corrections for masses, then the decay 
width of a Z boson into a neutrino pair is enhanced by the factor, 

d ) l-ji=A ^» 1 , / 1 a \ d I l-Ji=A ''i (Ar)X 

where Mz is the mass of the Z boson, and m is a neutrino's mass. If {1^, ■ ■ ■ , In+s) = 
0, then the factor ()4-2p certainly becomes unity. Dividing the summation in ()4-2p 
into two pieces, M| - 4m2 > and M| - 4m2 < 4 (f )' ^^+4* 

then the former part of the summation gives a positive finite value which is larger 
than unity, and the latter part of the summation becomes zero because the terms 
li = H and li = —n are canceled by considering Eq. (|A-3p . Thus, neutrinos cannot 
be thought of the fermions which cause the Casimir effect because the decay width 
of the Z boson is terribly enhanced. 

§5. Conclusions 

It has been investigated whether or not the Casimir effect from the fermions and 
the graviton can explain dark energy if they are the only massless particles which 
can go through the compact extra dimensions. In the second section, the finite tem- 
perature Casimir effect has been calculated on x (S"^)" space. And it has been 
seen that zero temperature Casimir energy and finite temperature Casimir energy 
can explain dark energy and dark matter, respectively, by considering the FL equa- 
tions on X (S^)^ X R space-time in the third section. The behavior of the energy 
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come from the Casimir effect may resemble that of tracker solution in quintessence 
modelP^ However, neutrinos cannot be the fermions which cause accelerating ex- 
pansion of the Universe, because the effect of extra dimension enhances the decay 
width of the Z boson. 

The forth coming problem is understanding the time evolution of scale factors in 
the era, which corresponding to matter dominant era in ylCDM model by assuming 
realistic initial conditions concretely. And it should be investigated that the effect 
to the matter density perturbation. 
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Appendix 



Aj{s,d,f3) defined in Eq. (|2-12p can be arranged by using the Abel-Plana for- 
mula, 

poo -I poo Jj. 

Y^F{n)- dtF{t) = -F(0) + i / [F{^t) - F(-zt)] , (A-l) 



n=0 

as follows. 



Z\j(s,d,/3) =d^ 



(27r)3+i 



(Z5+j,--W3)eZ"-^-i \ 
dt 



X sm 



\ 



where the formula, 
'A 



of {it) - Vit) = 2fe"^'^(*'-^') 



27r\ 
T J 



2 n+3 



sm['Ka]0{t - A), 



G7(z) = exp[aln(^"+z^)], 



(A-2) 



(A-3) 



has been used for A = 4^Jm'^ + + (27r/(i)2 ^"+5+^- tj. The factor l/(e^ - 1) 



I Z^i=5+j ''i 

expanded with respect to e~^ as follows to integrate with respect to k and t: 

d^+^k 



IS 



Z\j(s,d,/3) =d^ 



(27r)3+i 



E 



2d 



(Zs+j,--- ,ln + 3)(. 



TT \ 

yn-i-\ ^ 



r, \ 2 n+3 



i=5+i 
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u=l 



X sm 



s/3 



2 n.+3 



\d J 



i=5+j 



(A-4) 



Furthermore, the expansion of sinusoidal function with respect to its argument en- 
ables us to examine t integration: 



Aj{s,d,f3) =-d^ 



d^+^k 



(27r)3+i 



E EE 



n-j-l u=l v=l 



2d 



2s2/32y r(i> + i/2) J 2 



7r3/2s/3 \ ud ) {2v - 1)! 



X 2 n+3 1 2 

+ ^= + (1 E? 



i=5+i 



X K„ 



ud 



m 



^ + f E? 



(A-5) 



where K^{z) is the function called the modified Bessel function or Macdonald func- 
tion defined in Eq. (j2-6p . By transforming the Cartesian coordinate into the spherical 
coordinate in k integration gives, 



(«5+j,-A+3)eZ"-i-i«=i^'=i 



2s2/32y r(z; + l/2) J 2 /^2^ 



2d A/"^^ + (2vr/d)^Er=5+,^? 



3+2 
2 



2iTud 



ud J {2v-l)l ] '^['d 



2 n+3 • 2 

2 



Ee 



/27r 



2 n+3 



Mt) E'I 



i=5+j 



2(i 



Yl ■^3/2^^ 

a5+i,-A+3)eZ"-j-i«=i 



2 + (2^/d)2 Z2 



3+2 
2 



2'iTud 



v=l I ^ ^ 4=5+ 



2 n+3 > 2"^ 4 



5+i 



^=,2,2{^2_,(^)2j^n+3^^,,2| 
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Yl ^3/2^^ 

7n-j-l u=l ' 



3+1 



2TTud 



i+3 



X 7r2 (ud)" 2^ (m^ + [ 



2 n+3 
i=5+j 



d 



5 Ins 



2 



2 + (27r/d)2 Er=t+, ifV^WT^ 



i+3 

2 



(Z5+j,---,'n+3)eZ"-J-lw=l 



m2 + (2vr/d)2 Y.tL, n 



5+2 



27r 



J+5 
2 



m2 + (2^/d)2 Y.l=Lj l^iV^W+^ 



j+5 
2 



(A-6) 
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